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ABSTRACT 


In this thesis we report the results of a study of 
fitting of repulsive potential -energy (PF) curves for a 
few diatomic species and the ground and the first excited 
state surfaces for and of a qua siclassical trajectory 
(OCT) study of a model six -center exchange reaction. 

We have made a comparative study of various analytic 

y3 + 

functions for fitting repulsive curves, using H 2 ( 2^) as a 

test case. Exponential -poly function is found to be the 

best. Although the anti-Morse function is not as good as 

the exponential -poly, the former offers an advantage over 

the latter in that it could be extended to higher dimensions. 

We have found the rotated-anti-Morse-curve-spline to give a 

good representation of the potential -energy surface (PES) 

* 

for for linear as well as fixed nonlinear geometries. 

We have tested the utility of the Sorbie-Murrell (SM) 
function to fit the ground and the first excited state PES 
of using the chemically accurate PE values of Slegbahn 
and Liu for the ground state and slightly less accurate 
values of Nager and JUngen for the excited state. Static 
properties like minimum-energy-path, classical barrier 
height, force constant etc. of the ground state SM surface 
are in good agreement with those of the Truhlar -Horowitz (SLTH) 
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fit. Results of time-dependent quantum mechanical scatter- 
ing calculations for the collinear geometry using the 
surface are also in accord with the results obtained on the 
SLTH surface. 

We report, for the first time, results of a quasi- 
classical trajectory study of a model six-center exchange 
reaction 


over a wide range of relative translational energies {E,^) 

and for different vibrational states (v) . We find that 

, - 1 , 

at low (90-130 kcal.mol ), a four-center exchange 
reaction takes place with one of the molecules of the dimer 
acting as a mediating partner. At relatively higher 
(135-155 kcal*mol) ,for v = 8-11 of the dimer and the non- 
dimer molecules, a 'proper' six-center exchange reaction, 
with three/two/one stable molecule (s) as the products is 
observed. On the whole, we find that v is more effective 
than Eji m causing the reaction. 
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CHAPTER ONE 


INTRODUCTION 

1 -4 

molecular dynamical calculation is carried out 
by solving quantal or classical equations of motion for 
the nuclei on an adiabatic potent lal -energy surface (PES) 
and hence requires a detailed knowledge of the interaction 
potential. This is obtained by solving the electronic 
Schrodinger equation, within the Born-Oppenheimer approxi- 
mation, for a series of nuclear geometries, on an electronic 
computer. Although this can be achieved almost routinely 

5 

by using some of the readily available software packages 
like Gaussian 80 for a variety of systems, the enormity of 
the task of getting an ab initio PES of chemical accuracy 
(+ 1 kcal.mol can be realised from the fact that the 
number of systems for which such surfaces are available 
(for example, see ref. 6-9) can be counted by the fingers 
in one hand. Further, results of such calculations are 
available in the form of tables of numbers and for any 
subsequent dynamical calculation, these numbers have to be 
interpolated (and extrapolated if necessary) so that the 
potential and its derivative (S) may be computed at any 
arbitrary geometry during the simulation of the molecular 
collision. Here lies a 'bottle neck' in going from an 
ab initio electronic structure calculation to a molecular 
reaction dynamical study. 
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To fit the potential energy (PE) values in one 

dimension (ID - one independent variable^ for instance, 

bond distance for a diatom) is fairly straightforward. The 

variety of analytical functions and numerical methods that 

are available in the literature have been reviewed else- 

where^^'^^. Yet, one finds that fitting a PE curve with no 

minimum in it can be frustrating. Therefore, we have 
12 

undertaken a comparative study of various different 

*3 + 13 

repulsive analytic functions using H { E ) as a test 

2 

case. We have also tried the different functional fits for 
the ab initio potential data for Fe 2 ^^/ HeH^^ and NaLi^^. 
Details of the different functions used and a summary of our 
findings are given in Chapter 2. 

Pitting PE values as a function of two Independent 

variables (two dimensions ; 2D) is more difficult than fitting 

a function to data in ID. The choice of the functional form 

to be used for fitting is decided by the nature of the problem 

and the coordinates in which the data are available. The PE 

values for coll inear and fixed nonlinear geometries were 
17 * 

available for (2pE') in rotated-Morsp-curve-spline (RMCS) 
coordinates. Therefore we have tested the utility of a 
rotated-anti-Morse-curve-spline (RAMCS) approach to fit the 
purely repulsive potential (2pE' ) . Details of the 
fitting procedure and the parameters of the resulting fits 
are included in Chapter 2. 
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Fitting data in three dimensions (3D) is a niuch more 

difficult task and the state of the art is presented else- 

where^*^'^^. We only wish to point out here that the 

18 

Sorbie-Murrell (SM) function has been found to be successful 

1 Q 70 7 1 

in fitting ab initio PES for several systems: , H 2 O ' , 

22 22 +23 +24 +25 

LiPH , FH 2 , NH 2 , ChJ ’ 

Interestingly, its utility in fitting the most accurate 

8 

ab initio PES of Siegbahn and Liu^ for has not been 

2 6 

tested. While it is true that Truhlar and Horowitz have 

provided an excellent analytical representation of the surface 

and that its quality of fitting, with a root-mean-square (rms) 

deviation of 0.17 kcal.mol"^ may remain unsurpassed for some 

time to come, it would be admitted that it is not general in 

nature. Therefore, we have undertaken to test the utility of 

the SM function in fitting the Siegbahn-Liu PES. We have 

also tried to use it in fitting the PES for the lowest 

excited state (2pE') of by using the ab initio data of 

27 

Nager and Jungen . We have scrutinised the quality of the 
ground state fit by examining static properties like, barrier 
height, saddle point location, minimum energy path etc and 
also by comparing the results of a time- dependent quantum 
mechanical (TDOM) calculation of the probability for the 
Gollinear reaction 

H + H 2 -»• H 2 + H (Rl) 

on the SM surface with those on the Slegbahn-Liu-Truhlar- 
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Horowitz (SLTH) PES. The details are given in Chapter 3. 

The original aim of the present study was also to 

19 

use the approach of Varandas and Murrell via the QM 
function to obtain a double-valued functional representation 
of the PES for the ground and the first excited state of 

Some aspects of this attempt are included in Chapter 3. 
Traditionally, reactions like 

H 2 + I 2 '»■ 2HI (R2) 

were thought to go via a four -center (4C) transition state. 

But this picture started changing in the 60 's. The above 

28 

reaction was shown to proceed via an atomic mechanism. 

For the hydrogen -deuterium exchange reaction 

H 2 + D 2 2HD (R3) 

29-32 

shock-tube studies yielded an activation energy (E^^) 

of 42 ]<nal .mol~^ . The results were explained by a 4C 

mechanism postulating a vibrational excitation of the reactant 

33 

molecules, but the recent laser-beam experimental study 
of the reverse reaction 

HD + HD + ^2 

has shown that this 4C reaction does not occur even when 
one of the reactant molecules is excited up to the vibrational 


state v = 5. 



Extensive ab initio calculations (for example, 

34 

see ref. and references therein) yielded a classical 

barrier heighb (E^) £or the 4C reaction that was larger 

than the bond dissociation energy of H 2 . Detailed quasi- 

classical trajectory (OCT) calculations of Brovm and 
35 

Silver also failed to produce any exchange products 
via the 4C pathway. 

3 6 

Vtood ward -Hoffman rules also suggested that the 

reactions (R2) and (R3) may not proceed through a 4C tran- 
sition state. For such a pathway, the electronic ground 
state of the reactants is not correlated with the 

electronic ground state of the products and as a result, 

37 38 

a large barrier for the reaction may result. Wright ' 

suggested a six-center (6C) mechanism as an alternative 

as his LCAO-MO-SCF calculations yielded an of 90 heal. 
-1 

mol, significantly below the dissociation energy of 

39 

Thompson and Suzukawa used a semiempirical approach and 

-1 

obtained an Ej^ of 88 kcal.mol for the 6C pathway. 

, 40 

Dixon etal used a large basis set for their ab initio 
calculation and obtained an of 69 kcal.mol”^. With a 
still larger basis set, the of the reaction may get 
lower still I 

Crossed molecular beam experiments of Dixon and 
41 42 

Herschbach ' gave further support to the 6C mechanism. 



The reaction 


CI 2 + Br 2 2 BrCl (R5) 

did not occur at collision energies (F^) around 

Or -1 

^5 kcal.mol .However, at high pressures and low temper- 
atures in the supersonic nozzle expansion, where there was 
a high concentration of ( 012)2 they found the reaction to 
proceed at 1!]^ as low as 3 kcal.mol"’^ They attributed 
their results to 

CI 2 Cl 2 +Br 2 '>■ CI 2 + 2 BrCl 

BrCl....Cl 2 + BrCl 

(R6) 

involving a 6C transition state. 

We have undertaken a detailed OCT study of a model 

6C exchange reaction. We have used the diatomics-in- 

molecules approach to construct the PES- The parameters 

were adjusted so that Ej^ has the same value as reported 

40 

by Dixon et al . We have computed the reaction cross 

R 1 

section {0 ) over a wide range of E^: 90-135 kcal.T©!” 

for different v states of the diatoms. Details of the 
method SjUd the results are given in Chapter 4, 

A summary of our work and the conclusions are 
presented in Chapter 5. 



CHAPTER TWO 


PITTING REPULSIVE POTENTIAL -ENERGY CURVES AND SURFACES 

2.1 Repulsive Curves 

The first attempt to represent a repulsive PE curve 
was perhaps by Sato^^ and was for the triplet state (^2^) 
of H 2 . He had used an anti-Morse function of the form 

V = ^{exp(-2 3p) + 2 exp (-pp)},p = r - r^ (1) 

where D, r and f? are the bond dissociation energy, 
equilibrium bond distance and the curvature parameter 
respectively for the ground state ^2 have no 

such physical meaning for the triplet state. 

44 

Pedersen and Porter had extended the function by 
matching it with an exponential function at a longer range 

V = A {e'^fiP + 2e"*^^},r <1.6 au 

— rvF 

= b (r+c) e ,r >1.6 au (2) 

and were able to fit the ab initio potentials for H 2 and 

^2 Kuntz and Roach^*^ modified the anti-Morse 

function further by allowing f* to be a variable! 

f3 3 (1 typ + 6p^) (3) 

+ 45 48 

and fitted the potentials for H 2 and H 2 successfully. ' 

49 

Berces has shown that 



V 


(4) 


2 a 

= (D/a) exp (-2pp) {1 + (3 p) } 

is also more flexible than the anti-Morse in fitting the 

50 

triplet PE curve of ^ 2 , HI and I 2 . Kafri and Berry used 

an exponential "with a polynomial exponent for fitting the 

4 5 

ab nnitio data of Kolos and Wolniewicz 

V = exp(-(aQ + a^r + ^ 2 ^^ ^ a^r^)) 

(5) 

51 52 

Stine and Muckerman used ID spline for fitting the PE 

2 ■}* 53 

curve for H 2 ( • Sonneleitner and Beckel have shown 

that rational fractions can be used to reproduce the poten- 
, 54 

tial Eor H 2 (2p0^) 

~ 55 

More recently Varandas and Brandao have proposed 
the use of an exponential multiplied by an r~^ 

V = A r”^ exp (-bp ) (6) 

3 

for alkali dimers in their S state. Almost a decade ago 

56 

Murrell and Sorbie had proposed the use of an extended 
Rydberg (ER) function 

2 

V = ~D exp (-aj^p)fl+aj^p + a 2 p + ....) (7) 

for f-ittinq potentials with a single minimum. But there is 

no reason why it cannot be used to fit a repulsive curve in 

57 

the same way as an anti -Morse function, Huxley et al have 
reported recently that an appended ER function 
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2 3 

V = A (1 + a^r + 3 - 2 ^ + a 3 r ) exp (-a^r ) -tanh (r-r^/2) 

, -6 -8 -lOv 

x(Cgr -t Cgr + C^Qr ) 

(8) 

performs better than the ER in fitting a repulsive potential 

with a van der Waals minimum for a large number of diatomic 

species. In principle, given adequate data, the PE values can 

be fitted by a polynomial or better still by a ID spline or 
5 8 

Akima interpolation. In practice, one prefers a simple 
analytic function with a few parameters to be determined. 

To the best of our knowledge, no comparative study 

iTas been made of the validity of these various functions in 

fitting the data for a single system. The obvious choice for 

3 + 

a test case would bo the 2 state of for which accurate 

45 

ab initio values are available . We have undertaken such a 

study and the results are summarized below. We must mention 

that we have only considered global functions in this study 

and not considered the use of piecewise interpolation methods, 

43 

We have specifically considered only the anti-Morse , 

49 50 55 

Berces , exponent ial~poly, modified exponential and 

the functions as we are interested in functions which 

require only limited input data and also are extendable to 

higher dimensions. 

45 3 + 

Ab initio values of InV for ( 2^) are plotted 
against r in fig. 1 and for comparison, the different 
analytically fitted potentials are shown in the same graph. 



V (eV) 



Fig.1 Comparison of different analytic fits to the ab 
initio data for H 2 fey). The ab initio PE 
values are connected by a solid line. The 
dashed lines represent the different tunc- 
-tional fits • 




The parameters for each function were determined by a 

nonlinear least squares fit of the ab initio data using the 

59 

B04FCF program of NAG library . The rms and the maximum 
deviation for each function are listed in Table 1. 

3he anti-Morse function fits the jpotential adequately 
at short range. The modification introduced by Berces 
improves the fit dramatically - but at the expense of adding 
a parameter. The exponent lal -poly function deviates much 
more than the Berces function at long range. But the former 
gives the best overall fit as is illustrated with the aid of 
a residual plot in fig. 2. The ER functional fit is only 
slightly inferior to the Berces, but involves as many para- 
meters as the exponential -poly fit. The modified exponential 
function gives the poorest fit. 

We have tried the different functional fits for the 
ab initio potentials for He 2 y HeH and NaLi.^^”^^ The results 
summarized in Table 1 show that while the anti-Morse fuction 
fits the potential reasonably, the exponential-poly gives a 
better fit for all the systems. The quality of the fit 
obtained by using the extended Rydberg function is comparable 
to that of the exponential -poly for He 2 and NaLi. Berces 
gives the worst fit for NaLi and extended Rydberg that for 
HeH, If a large number of potential -energy values are avail- 
able and the number of parameters can be as large as six, the 
exponent lal -poly fit is obviously the best or close to being 
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the best. 

It also offers the advantage that InV can be fitted 
to a polynomial by a linear -least -squares fit to yield an 
excellent set of initial guess parameters for carrying out 
a nonlinear-least-squares fit-an important practical 
advantage over the other functional approaches. On the other 
hand, if the number of parameters has to be as low as th^ree, 
the anti-Morse function yields reasonable results and it also 
promises to be the candidate when PI? values involving more 
than one independent variable are to be fitted. 

2.2 Repulsive Surfaces 

As a test case of a repulsive PES in 2-3 independent 
variables, we have examined the possibility of fitting the 

if 

PES for (2pE') in its geometry . This surface is . 

needed in the study of the spectroscopy of the transition I 

I 

state in predicting the wings to the Lyman-a line in 

6 0 I 

H + collisions . The amount of ab initio information ' 

61 I 

available on this surface is scanty . Therefore, at our ! 

I 

1 7 ’ 

request, Raynor has computed a diatomic -in-molecules PES 

6 2 I 

in rotated Morse coordinates (A*0) (illustrated in fig. 3) 

j 

for linear as well as nonlinear geometries i three PE values | 

I 

along each ray, for 0=0, 10, 20, 30, 35, 40, 42, 43, 44, j 

I 

45^^ and hence by symmetry at O = 46, 47, 48, 50, 55, 60, 70, | 

80, 90° for the bending angle y = 60/ 00, 120 and 180° . i 


1 . 




io 


The idea of using a rotated Morse coordinate is that 

for 0- = 0°, we have the PE curve for the reactant diatom 

and by rotating it about the swing point S(R°, R^) through 

0- yields the PES, the 0 = 90° limit providing the PE curve 

for the product diatom. In principle, the Morse curve 

representation can be replaced by any other suitable repre- 

* 

sentation. For only three data points are available 

along each ray and therefore, we were constrained to use 
only 3 -parameter functions in our rotated-f unction approach. 
We tried to use 


and 


2 

V ~ exp (a^ -I a^p + a 2 P ) 

V = D£ exp 


(9) 

( 10 ) 


which are the truncated and modified forms of the exponential- 
poly and ER functions respectively. We did not obtain 
satisfactory fits at low energies using eq. (9) . Although 
eq. (10) gave a good fit along each ray, the variation of the 
parameters (D, a^, 32 ) with 0 was not smooth. Pence we were 
not in a position to fit the potential successfully using 
these two functions. 

★ 

For the col linear geometry of the 0 dependence of 

D,f! and Jl_ are shown in fig. 4. These were subsequently 
fitted by ID spline interpolation to obtain the overall fit. 



D(eV) P(A ) le(A) 



0 (deg) 

Fig. 4 Dependence on 0 of the anti- Morse 
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The resulting PES is illustrated with the aid of PE contours 
in fig. 5, We were also able to obtain a similar fit to the 
potential for bent geometries The 0-dependence of the 

parameters for y - 120, 90, 60° is illustrated in figs. 6-8 
and It is comparable to that of Y = 180° except for a double 
hump in for Y = 60°. The overall fit for Y = 120, 90, 

60° IS illustrated with the aid of PE contours in figs. 9-11. 
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Fig. 9 Same as fig* 5 for Y" = 120°. 
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Fig. 10 Same as fig. 5 for Y" =90°. 
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CPAPTER THREE 


SORB IE -MURRELL FUNCTIONAL FIT FOR THE PES FOR THE GROUND 
AND THE FIRST EXCITED STATE OF 

3 . 1 Ground State PFS 

For a triatonic system ABC, the SM functaonal repre- 
sentation of the PFS IS given by 

'^AEC ' V''b+V''aB»1> ''cA^s' + VjlR^.Rj.Rj) 

( 11 ) 

where Vg, are the monoatomic terms which can be set to 

zero, '^BC' ^CA 2 -body terms as a function of the 

bond distances Rj, / R 2 -- R 3 respectively and V^. is the 3 -body 

term that goes to zero as one of the {R ^ lvalues goes to 

infinity. We have used the FR function (eq. (7) in Chapter 2) 

to fit the 2-body potential, that is, the PE values of Kolos 
13 

and Wblniewicz for H 2 . The resulting parameters as obtained 

by using the nonlinear least-squares fitting subroutine 

59 

ED4fCF of the NAG library are reported in Ibble 2. 

The 3 -body interaction term is represented by 

= (P + P')TA (12) 

where P and P' are polynomials and TA is a range factor. 
Because of the symmetry of the system ^ use symmetry- 
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Table 2 • Paranietersfor the Extended Rydberg Functional 

£it for 


D /eV 

.o 

33/2-^ 

a4/X'^ 


4 747733 
0 74139 
3.6720 
2.9869 
2.7045 


-1.0972 



adapted displacement coordinates 




/l73 

/ i 73 

/I73 




^2 


0 

/T72 

-A/2 


^2 


O 3 


/273 

-/1/6 

-/I76 


^3 

(13) 


where Pj^ ~ ^i"^i' ^i reference point, treated as an 

adjustable parameter We use a fourth order polynomial 
for P : 

P = + a^Oi + or^Oi + 

®8^1^r ^ (14) 

where { " 0/ 10 at-e adjustable parameters and 

Pi 2” 

i (15) 

and 

^s “ 302^^3 (^0) 

1 9 

The additional polynomial P* was introduced to reproduce 
the Jahn-Teller intersection for the geometries, as 

P' = qj. {Pi + 320 ;l ‘' 3^1 **■ ^ 

(17) 

with 1 = 1-7, and the positive root of eq. (15), We 

took TA to be (l-tanh(Y0]^/2) ) m which y is yet another 
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parameter to be determined. 

Omission of the polynomial P' led to an rms deviation 

of 3.0 kcal.mol"^, and a maximum deviation of 15 kcal.mol~^ 

for noncollinear geometries. Therefore we included P* and 

obtained the 20 parameters by a nonlinear least squares fit 

of the function to the 267 PE values of Sieybahn and Liu^~^ 

The result was an analytic fit with an rms deviation of 
-1 

0.78 kcal .mol . But there was a spurious potential well of 

about 4 kcal .mol depi h at Rj - 2 . 1 2 5^ and R,, = 0.7 4 R in 

the entrance channel and a similar one in the exit channel 

for the colllnear conC Lcfuration Thjs was mainly due to 

inadequate ab initio data at large separations of H from H 2 - 

An addition of 14 Pg values, generated from the SLTF surface 

for geometries with = 2.12-2.91 R and = 0.71-0.77 R 

to the ab initio data resulted in a lowering of the well depth 

to 1 kcal,mol~^. Therefore we decided to add yet another 28 PE 

values generated from the SLTH surface in the same (Rj,R 2 ) 

range. We also added the exact location and energy of the 

saddle point of the SLTH surface to the data. Thus the function 

was fitted to a total of 310 input points and the resulting rms 

deviation was 0.825 kcal.mol"^ and the well depth was only 

0.50 kcal. mol The clu&sicol bcirraei" height, however, was 

-1 

only 9.42 kcal. mol . Therefore we added a correction term 
0.013 (I -tanh (0^ -J . 33 774 ) ) to the fund ion, to make the barrier 
height 9.8 kcal. mol We also added a gaussian in the form of 
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%uss " {-1.92 (0^ - 2.87796)^} (18) 

to reinove the well. The overall fit now has the correct 

—1 

classical barrier height oE 9 8 kcal.mol and an rms 

deviation of 0.9 kcal <nol . The well in the entrance channel 

1 

for the collinear configuration is only 0.12 kcal.niol" deep. 

For an approach angle of between H and 1*2, there is a 

well of 1.7 kcal.mol But it is unlikely to affect the 

dynamics as the reaction is dominated by the collinear geometry. 
The parameters oE the final fit are given in Table 3. 

Before proceeding to a detailed dynamical calculation it 
IS essential to ensure the quality of the fitted PES using 
criteria other than the rms deviation. ThereEore we have 
plotted the PE contours for collinear and different fixed non- 
collinear geometries in figs. 12-16. 

The variation of the classical barrier height for the 
reaction (Rl) with y is illustrated in fig. 17. For a 
comparafon, results obtained by the SLTil fat are Included in 
the same figure. The minimum energy path (MEP) and the 
variation oE the potential along the MEP as obtained by the SM 
fit, for the collinear conEormat ion, are illustrated in figs, 18 
and 19. We have computed the stretching force constants 
Kg = (3^E/3g^) and (a^E/Sg^) 'Cor collinear at the saddle 

point = g 2 = 0) where g^=(/3/2) (Rj^ + R 2 - 2R ) and 
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Table 3 • Poramelors for the RM f n L for the ground and the 

o 

first excited state of . Distances are in A 
and energy in eV units 


Parameter 


Ground State 


First Excited 
State 


«1 

a 2 

“3 


ag 

Otg 

“lo 

3 1 

32 

33 

34 

35 

36 

^7 

Y 


13.045 
1.7038 
- 0.79314 
- 0.89915 
0.89972 
- 9. 1007 
6.3079 
0.03498 
0.01479 
3.3866 
- 12.484 
- 6.6106 
- 2.3982 
2.7970 
13.547 
- 1.9411 
10.431 
- 7.6274 
1 .4739 
0.7937 


- 16.615 

- 8.2855 
- 15.927 
125.54 

- 2.7462 
- 2581.4 
- 57.961 
- 10.528 
704 28 

76. 219 
- 298.05 
90.794 
- 787.63 

- 5.4865 
- 265.62 

471,29 

- 1594.6 

- 78.548 

2.8683 

0.7937 


R2(A) 



Rl (X) 


Fig. 12 PE contours for the colli near (Y‘= 180°) confor- 
mation from the 5M fit. The numbers give 
the PE values in kcaLmoP^ units relative 
to the separated H and H2 - 




^•9 • 13 Same 





Rl (A) 

■fig • 12 for Y = 90°. 


Fig. 15 Same as 


R2(A) 



Fig. 16 Same as fig- 12 to 



Classical barrier height (k cal /mole) 



Fig. 17 Variation of the classical barrier height 
with T for SM and 5LTH surfaces- 



• » 



^Am\ 


for the collinear 
and SLTH fits- 


I 




Arc length (cm) 


Fig. 19 PE for the collinear reaction (R1) 
as a function of distance along 
the minimum-energy path . 
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02 = (1/2) (Rj^ - R 2 ) and compared with the values for the 
SLTH surface A summary of the saddle point properties is 
given an Table 4 . 

Prom the above figures it is clear that the fit does 
not have any unphysical behaviour except for the well of 
1 .7 heal. mol ^ at y - 90°. In addition it meets many other 
criteria required of an accurate analytic fit; 

( 1 ) It is symmetric with respect to an interchange of and R 2 - 
(ii) It has a correct barrier height for the collinear 

reaction and it fits the ab initio data satisfactorily. 
Asymptotically it gives an accurate description of the 
diatom . 

The angular variation of the potential as obtained by 

the SM fit IS illustrated in figs. 20 and 21 for r„ = 0,74 X 

and 1.48 A respectively. For comparison /the plot for 

r„ = 0.74 S obtained from the SLTH surface is shown in 
^2 

f ig . 22 . 

It is instructive to expand the potential as 

•V(Rj,R2.R3)= '’l 

where R 2 is one fixed bond distance and X 2 angle 
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Table 4 : Saddle Point Properties 


Property 

Siegbahn-Liu 

SLTH 

Present fit 


0.9297 

0.9297 

0.9432 

(kcal.mol"^)^ 

9 .80 

9.802 

9.80 

Kg (eV 

31.09 

31 .40 

34, 78 

, o_2 d 

Kg (eV A 

5.63 

-5.60 

-8.31 


a is the distance between adjacent hydrogen nuclei at 

the linear symmetric saddle point 

b classical barrier height taken with respect to 

the corresponding H and 112 results, 

c IS symmetric stretch constant. 


d is asymmetric stretch constant- 


2.4 



Fig. 20 


Fixed-molecule PE contours at 
for the 5M surface . 




2.4 



Fig. 21 


Same as fig. 20 for 



1.48 A . 




^15 


between the bond vector and the cent er -of -mass separation 
vector. At a given R„ „ the V. values are obtained by 

li n 2 JLj 

integrating V(r 5 ,r ) weighted by the appropriate 

^ 2 

Legendre polynomial 


(R- 






1 

/ 

-1 


Pt. (<^05X5) 


V(R2/Rjj_h^ 


)d (cosXo) 


( 20 ) 


We have computed and presented in Table 5 ,V^ and V 2 for 

Pu u = 1 • 6-3 .5 a^ . 

H-H 2 o 

Over the entire range, of the SM fit is in good agreement 

with that of the SLTH, of SM is also in agreement with 

that of SLTH at R„ „ = 1.6 a . But the former is larger in 

H-H 2 o 

magnitude than the latter at R„ „ = 3.5 a . Similar to the 

H-H 2 o 

SLTH fit, the SM fit also yields a V 2 that changes sign thrice 
along the R„ „ axis. It is positive for R„ „ =1.6 a_-2,la_ 

iTi — ^2 o O 

and becomes negative at S.la^ and becomes positive again at 

6 3 

2.9a^ similar to the results of Varandas . It becomes 

negative again at larger R„ „ . Ihe occurrence of the rniddle 

H-H 2 

zero of SLTH at a larger R„ „ need not be accurate and 
^ ri'-ri2 

ours Miay bo beLLor in Ihjs 1 ospocn 

The ultimate test of the reliability of a fitted surface 
would be to carry out dynamics on it and compare the results 
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Table 5 

: Legendre 

coefficients 


in kcal.mol 


units for 

the SM and the 

SLTH fits 




SM 


SLTH 

<®o> 






L=0 




1.6 


61 . 50 


61 .60 

2.0 


32 61 


31.67 

2.-5 


17.065 


16.12 

3.0 


8.867 


8 699 

3.5 


4.063 


4.477 

4.0 


1 290 


2.100 

4 . 5 


-0.124 


0.838 

5.0 

L=2 

-0 600 


0.266 

1.6 


25.92 


26,48 

1.8 


7. 473 


- 

2.0 


0.546 


-0.922 

2.1 


-0.831 



2.5 


-1.255 


-2.597 

2.8 


-0.102 



2.9 


0.291 


- 

3.0 


0.663 


-0.911 

3.5 


2.025 


-0.170 

4.0 


2.4 98 


-0.041 

4.5 


2.239 


-0.004 

5.0 


1.644 


-0.006 

6.4 


0.121 


-ve 

7.1 


-0.207 


- 

10.0 


-0.193 


-ve 
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with those on the exact surface Since a TDCM study of a 
collinear reaction can be carried out almost routinely 
(for example, see ref. {(60c)}) and the reaction (R1 ) is 
dominated by the collinear approach, we carried out a TDOM 
calculation on the EM surface and also on the SLTH surface 
for the collinear geometry. A comparison of the results 
on the two surfaces is presented in fig. 23 in the form of 
the reaction probability as a function of time for H 2 in 
Its first excited vibrational state at a translational tem- 
perature of 1000 K. It IS clear that the SM fit not only 

has the static properties similar to those of SLTH but 

\ 

It IS also able to reproduce the dynamical results as obtained 
on the SLTH. Therefore we can say that the SM fit is an 
adequate representation of the Siegbahn-Liu surface for the 
ground state of . 

3.2 Excited State PES 

A correlation diagram of the first few states of H^ 

as given in ref. 60(b) is reproduced in fig, 24. We have 

undertaken to Cit the available ab initio data Cor Lha 

lowest excited state (2pE' ) that conically intersects with 

27 

the ground state in the geometry. Nager and 0 ungen 
reported the binding energy (e) values for the third electron, 
in the frozen core approximation for the meta stable Rydberg 

■k 

states of upto the principal quantum number n = 4 in 
the neighborhood of the equilibrium geometry of the parent 



TIME STEPS 


Fig. 23 Reaction probability as a function of 
time with each time step equal to 
0.53875 X on 5LTH and 5M surfaces, 






Fig . 24 Correlation energy diagram for the four lowest potential- 
energy surfaces of H3 together with the ground -state 
surface , as a function of nuclear geometry • The 
lowest excited-state and ground-state surfaces meet 
in a conical intersection at 031-^ geometry . 



ion . For n =2, they reported values of the total 

-f 

energy {that is a sum of e and the core energy for H^) for 
Cg, C 2 ^ and geometries. For and geometries, 

however, they reported only the e values and we added the 

+ 64 

core energy for as obtained from the Giese-Gentry 

fit of the ab initio results of Csizmadia et al in order 

to obtain the total energy. In toto, there are 67 PE values 

available for the (i 2 ( 2 pE'). 

We followed essentially the same strategy as we 

employed for the ground state in fitting the SM function to 

the PE values Cor the tirst excited state of . Asympl'o- 

tically, H2{2pE‘) goes to H(ls) + used the 

exponent la i-poly function (eq (5) in Chapter 2) To fit 

1 3 

the ab initio PE values of Kolos and Wolniewicz for the 
repulsive state of n* ^ud the parameters are listed 

in Table 6. 

The 3 -body interaction potential was represented by 
eq.l2 with the P' defined in terms of the negative root for 
in eq. (15) . The resulting functional fit had an rms 
deviation of 3.1 kcal.mol'’^. An examination of the 
variation of the PE as a function of the side length of the 
equilateral triangle for °3h geometries as illustrated in f ig . 25' 
reveals the difference between the fits for the ground and 
the excited state PES of . The two fits do not yield 
identical PE values for the 03h geometry. Part of the reason 



i; ' I-'*/'... ' 

/ / 1,1^ nn }r 


»7m* 


. «. A ^9209 


r,l 


Tci bl e 6 


Parameter b Cor the exponential-poly Lit 

*■ 3 H- 

of the respulsive PE curve for 


^o 

4.9986 

ai/ 8 "^ 

-8.3407 


9.0983 

a 3 /A 

-6.4784 

,0-4 

34 /^ 

2.1694 

,0-5 
a 5 / A 

-0.2802 



Fig. 25 PE in kcal,mor’ units as a function of the side 
length of the 03 ^ geometry R(X) for the 5M 
fitted ground and excited state PESs . 







for the discrepancy is the difference in the quality of 

6 —8 

the results of Siegbahn and L3u° results and those of 

27 

Nager and Jungen . Since the former values are more 
accurate than the latter, we replaced the Nager-Jungen PE 
values for twelve geometries by the values generated 

from our SM fit of the Siegbahn-Liu surface. The resulting 
fit had an acceptable rms deviation of 2 6 heal. mol The 
parameters for the fit are included in Table 3, 

3.3 Attempts to fit a Double-Valued Surface 

19 

We have tried the approach of Varandas and Murrell 
of using the SM function as an analytic fit for both surfaces 
as two roots of a polynomial. The parameters a^, a^, 

a^, arj and y were obtained by fitting seven ab initio values 

7 

of Siegbahn and Liu . We attempted to obtain the remaining 
twelve parameters by fitting the function to the data of 
Nager and Jungen for and geometries (56 points in 

all). But the parameters did not converge. 

For systems like 0H2^^ and the SM 

approach has been extended to fit the double valued surfaces. 

The basic philosophy has been to fit each diabat to an SM 
function and the nonadiabatic interaction term to a trigono- 
metric function and to obtain the lower and upper surfaces 
as the roots of a 2x2 determinant. Such approach would be 
successful if equal number of ab initio values are available 



£or both the surfaces. Unfortunately, sufficient number 
of data points are not available for the excited PES at 
this stage. V/or'k is in progress to arrive at a double valued 
function for 



CHAPTER FOUR 


DYNAMICS OF A MODEL SIX-CENTER EXCHANGE REACTION 


A quasiclassical trajectory study of a Tiodel six-center 
exchange reaction between a van dec Waals (VDW) dimer 
A 2 . . . B 2 and a monomer 
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C 

I 

c 


^A 

A^ 


B 




'B' 


/ 


AB J AC + BC (R7) 


was initiated earlier in our group (The study yielded only 
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nonreactive trajectories) . We have followed it up to 
obtain results like reaction cross section and rate constant 
data which can be compared with the available experimental 
results and also an insight into the details of the dynamical 
pathway. We give below a brief description of the methodology. 


4.1 Methodology 

We have used three types of coordinate systems (each 
being cartesian) in our study. 

( 1 ) Space-fixed coordinates XYZ with the center -of -mass 

(c.m.) of the whole hexa -atomic system as the origin. 



(ii) Space-fixed coordinates X'Y'Z' with the c.ta. of the 
tetra-atomic species (VDW molecule) as the origin. 

(ill) Molecule-faxed coordinates xyz for A 2 / x'y'z' for B 2 
and x"y"z" for C 2 with the c.m. of the respective 
molecule as the origin. 

Figure 26 shows the coordanates and their conjugate 
momenta in the space-fixed frame XYZ and Fig. 27 the space-fixed 
X'Y'Z'. 

We assume that the motion of the atoms can be described 
70 

by Hamilton’s equations. In order to solve these equations 
we need the derivatives of the potential with respect to the 
coordinates in addition to the initial conditions (coordinates 
and momenta ) . 

71 72 

Based on the London equation, Eyring and coworkers ' 
derived the secular equation for a six-electron system to 
evaluate the PE : 

- E = 0 (21) 

where = ‘^Xj^lx^^ with fhe basis 

■^9 

functions. Thompson and Suzukawa applied the same approach 

to the 6C reaction. We have extended it to get a value of 

40 

equal to that of the ab initio value of Dixon et al 

We have evaluated the individual terms in the secular 
7 1 

determinant which consists of Coulombic integral and 
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Fig . 26 Coordinates and their conjugate 

momenta in the space-fixed frame XYZ. 



2 ' 



Fig. 27 Orientation of the VDW bond and 
the coordinates of the c.m. of A2 
and 02 in the frame . 



exchange integral J^j/ which are functions of the corres- 
ponding internuclear distances. We have used the London- 

73 

Eynng-Polanyi-Sato (LEPS) formalism to evaluate these 
integrals 


1 - 


1 3 

where E and E are the sing] et and triplet energies for 


the diatomic molecule given by 


Eij = lexp{-2a^j - r^^)} 

-2 exp (r^j - (24) 

e 

IS the equilibrium internuclear distance and is the 

Morse parameter. The triplet state curve for the atom pair 

. 44 

ij is represented by the Pedersen -Porter function 

^®ij (exp{-2fi^j (r^j-rJ^ ) } + 2 j ) > ) 


for Rp and 




for r^j> 


( 26 ) 



OO 


7 4 

The parameters computed by Raff et al for H 2 
been used in a number of trajectory calculations for three- 
and four-atom systems results of which were in satisfactory 
agreement with the experi'^iental results Hence we made use 
of these parameters, except for the Sato parameter S, which 
was varied till for the reaction (R7) through a regular 
hexagonal geomotry agreed with the ab initio value of 
ty kca L .moi“^.The parameters are listed in Table 7. 


The derivative of E with respect to any specific 

internuclear distance has been calculated by making use of 

75 76 

the Hellmann -Feynman theorem (HPT) . It can be expressed 
a s 




3^1 



-> 





S Cj^ 


(27) 


where are the eigenvectors corresponding to the lowest 

eigenvalue Cj^ and Ej^ were calculated by using the 

77 

NICER subroutine. These derivatives were then converted 
into derivatives with respect to the coordinates using the 
chain rule. 


69 

Initial coordinates and momenta of the molecules 
involved in the reaction were calculated as follows. Initially 
the molecule A 2 was aligned along the Z-axis with its c.m. 
as the origin. Hence in the molecule-fixed system the 
coordinates of A 2 are given by 
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Table 7 i Potential-energy Parameters for 


1 -1 
D (kcal.mol ) 


109 46086 



0.74127 

1 97357 

45.35618 

1 88999 
1113.59090 
0.529167 
3.166560 
0,84667 
0.32000 
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Xi - - X2 = y2 - 0 

^2 “ ~^i/2 (28) 

where is the internuclear distance of A 2 - The initial 

x,y,z components of the total momentum conjugate to the 
molecule-fixed coordinates are 

= TiiVy<''2> 

'’y,^ = 

= |Pvl 

Px2 = ™t2“y<''2* 

Py^ = (A2) 

p = - j p I (30) 

=^2 ' vl 

where to (A^) and to (A-,) are the angular velocities of A^ m 
X and y directions respectively and is the vibrational 
momentum. 

The initial coordinates and momenta of the molecule in 
any position in 3D-space can be calculated by multiplying by 
the appropriate combination of the rotation matrices 
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R,(g) 


Ry(3) 


cos 3 
-sin 3 
0 

cos 3 

0 

sin p 


sinP 

cos 3 

0 

0 

1 

0 


0 

0 

1 


~sin3 

0 

cos 3 


( 3 i) 


The coordinates of the niolecule a.^{x,y,z) and of r (x' 

2 ^2 

transformed into space-fixed coordinates (X'Y'Z') can be 
evaluated by the following transformations (see fig. 27 } 


y ' z ' ) 



tr: 

P3 sinO-,^ COS 04 
pj sin 9 ^ srn^^ 

+ (0il.Ry(O^) 

^1 

^1 



Pjl cosG^ 


- 


— - 


— — 


r 



-P 2 sinO^ cos0^ 


K' 

1 


- 

-P2 sin&^ sin0^ 

■> Rj^'2)Ry(02) 

y ' 

Z ' 

X 


-p , cosO ^ 

* 2 A 

( 

2 ' 

X 


l:=l,2 


( 32 a) 

1 = 3 , 4 

( 32 b) 


021 and ^3 are the orientation angles of A2, B2 and 
C2 in the molecule-fixed frames xyz, x'y'z' and x"y"z" 
respective] y. P2 the distances between the c.m. of A2B2 

and that of A2 and the c.m. of A2B2 that of 63 respectively. 



^4 orientation angles for the VDW bond, 

of da stance p 

The corresponding momenta in the space-fixed frame 
wall bo given by 




p 

1 

X ’ 

1 



1 


Py 
^ 1 

P 

z 

1 


1 


p 

xi 



Py. 

1 



Pz' 


Pz' 

^ ' 


^1 


1 , 2 


3 , 4 


( 33 ) 


The collision between (A2B2) and C2 was studied by transfor- 
mang the coordinates X'Y'Z' to XYZ. If the c.m. of (A2B2C2) 
IS taken as the origin of the XYZ frame, we have 

2(mj+m2) c + ^ 

where m2 and m^ are the atomic masses of B and C respectively, 
V = v-c is assumed to be the relative velocity, as v is the 
c.m. velocity of C2 and c the c.m. velocity of (A2E2) , lies 
along the Z-direction. Hence it follows that 



bS 


kj+m^) 


C + 1113 (V + c) = 

^ -ni3 V 

n = — = 

^1-I’D2+'t»3 


and 


(nf)j+Tn2)V 
'll j ^tn2^ 1TI3 


0 


The position and momenta of the €3 molecule in the XYZ frame 
(Pig. 28 ) are; 




"‘'a 


,,((( 3 ) Ry„(»3) 


X. 


y 


JL 


Z " 
X 


P 


X 


II 

1 



6 


( 34 ) 


where b is the impact parameter and the distance between 
the c.m. of €3 and the c.m. of (A2B2). In the above trans- 
formation, Z IS assumed to be the direction of approach of 
the C2 molecule to the (A2B2) system and the c.m. of (A2B2) 
and C2 lie in the Y 2 plane. 
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In the xyz coordinate frame the initial coordinates 
and momenta are given as follows : 


— 


— — 



X 


0 


X' 

1 




X 



mib 



Y 

1 


(mj^+m 2 +m-^ ) 

4* 







1 


(mj^+m2+i^'3) 





0 



1 




1 



0 

+ 


Y 




Y' 

1 




1 



'V 3 V 



1 — 

•u 






mj^ “ when i = 1 and 2, and = '^2 1 = 3 and 4. 

In the present study we have fixed the values of the 

vibrational and rotational quantum numbers v, J and the 

relative translational energy (E^) while b has been selected 

nn a stralifJed manner from the appropriate dJsLribution 

function. The orientation angles have been selected randomly 

in the interval generated by the IMSL subroutine 

7 8 

GGUB . After selecting initial conditions corresponding to 
coordinates and momenta the equations of motion were solved 

using the fifth-order predictor -corrector subroutine adapted 

79 ”16 

from the REACTS package' . A step-size of t= (0 . 5-1 . 0) xlO s 



Wd& used for Lhe intecjraf .ion 


We analysed the final state of the trajectories as 
described below. In the collision between (A 2 B 2 ) and C 2 the 
possible channels are : 

')• A2 + B2 + C2 (R8) 

-5- AB + AC + BC (R9) 

-»• 2An + C 2 (RIO) 

2AC ■) B 2 (RIl) 

->■ 2BC + A 2 (R12) 

69 

(R8) is a nonreactive process, (R9) a 6C reaction and in 
(RIO), (Rll) and (R12), C 2 / B 2 and A 2 act as mediating partners 
respectively. The various possible products formed during 
the reaction, along with their identifications are listed in 
Table 8. 

To identify which products have been formed, we made use 
of the internuclear distance criterion. All the fifteen 
internuclear distances were compared and the integration was 
continued until only three of the fifteen distances were less 
than a specific distance R^. The molecules corresponding to 
these three distances were identified and their internal 
energies analysed. 




69 


bie 8 


Species Labels for Different Reactive Channels 


Process 


Coinbination of atoins in sets 
of three 


Nonreactive (R8 ) 

12, 

34, 

56 






AA, 

EB, 

CC 




Six-centered reaction 

13, 

25, 

46 

13, 

26, 

45 

(R9) 

AB, 

AC, 

IC 

AB, 

AC, 

EC 



14, 

25, 

36 

14, 

26, 

45 



AB, 

AC, 

EC 

AB, 

AC, 

BC 



15, 

23, 

46 

15, 

24, 

36 



AC, 

AB, 

bC 

AC, 

AB, 

BC 



16, 

23 , 

45 

16, 

24, 

35 



AC, 

AB, 

BC 

AC, 

AB, 

BC 

C 2 acts 

as a niedDating 

13, 

24, 

56 

14, 

23, 

56 

partner 

(lUO) 

AE, 

AB, 

CC 

AB, 

AB, 

CC 

B 2 acts 

as a fnediaLincj 

15, 

26, 

34 

16, 

25, 

34 

partner 

(Rll) 

AC, 

AC, 

BB 

PC, 

AC, 

BB 

J ^2 act s 

as a mediating 

12, 

35, 

46 

12, 

36, 

45 

partner 

(R12) 

AA, 

ry^ 

BC 

AA, 

B2, 

BC 
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The c m . velocity in the space-fixed frame for product 
molecules was calculated as follows 


(V 


ij 


'X 



+ ) 


(V 


IJ 


Y 




+ m^ ) 


1 




(p. 


>1 * 




(m^ + nij ) 


( 36 ) 


x<y/Z components of the linear momentum in the molecule-fixed 
coordinates are 


X, 


= P. 


X- 






= P 


Yk 




P' 

Z 1 




( 37 ) 


The internal angular momentum is calculated by transforming 
the space-fixed co-ordinates to molecule-f ixed coordinates 


as follows : 
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(ni^ + .n X ) 






^’^1 ^1 ■*' *^1 ^1 ^ 


(m^ + nij ) 


z, = Z, 
k k 


(m , Z , + m Z ) 


(td^ + m^ ) 


/ k = i,j 


For diatomac species ij , the final internal angular 
momentum is given by 


= (vi ) + ivj p; - “j p; > 

X J J 


(Mlj)y “ ki Px^ - Xi PjJ + (Zj Px^ - Xj Pz^> 


= ‘>‘1 Py, - Px > + Py - S'j Px ) 

*^1 X 11 


2 

The corresponding rotational energy E^Q_j.(ij) = /2I 

2 2 
where I = yr.. with vi= ni m / (m +m ) and r ={ (x -x ) + 

J-J XJ XJ XJ XJ 

^ and ‘"ij’y'* 

The vibrational energy of the product molecule (i1 ) is 


calculated Erom 



E' 
vib 


+ V,, - (13) 


xj rot 


(40) 


where is the internal kinetic energy: 


= <l/2) { (P/ P'^ + Pi^)/m 

1 1 




( 41 ) 


and the Morse potential. The final translational energy 




is 


^ ^TOT ^ib ^rot 

4.2 Results and Discussion 

For initialising the tra3ectories, we set the VDW 

, o o 

bond distance Pj^ + ~ 3,0 A and = 6.5 A. We set J = 0 

for all the trajectories We chose fixed values of v and 
as described below. Trajectories were computed for 
different b-straia with b chosen differentJy Cor 
different (v, Generally, about 100-200 trajectories 

were run in each stratum. On the whole, about 3000 trajectories 
were run, with each trajectory taking about 3 min cpu time on 
the DEC-1090 computer Details of the number of the traject- 
ories run tor different initial conditions and the reaction 



R 

cross section (a ) are given in Table 9 

The accuracy of integration of each trajectory was 

monitored by energy and angular momentum conservation 

measures. In many test cases, we plotted the different bond 

distances as a function of time to determine the outcome of 

the collision. This helped us in ascertaining the mechanism 

of the reaction as well. We find that for (v^ ,Vg ,v^ )=:(0,0,0) 

1 2 2 2 
at = 90-135 kcal.mol" , the reactions (Rll) and (R12) 

take place through a 4C pathway with the third diatom acting 
as a mediating partner as is illustrated for a typical traje- 
ctory Ln fig. 29, For higher vibrational states like (11,8,8) 
at = 13 0-150 kcal.mol~\ we observe a 6C reaction. That 
is, the three (reactant) diatoms lead to three different 
(product) diatoms as is shown in fig 30 for a single tiaje- 
ctory. At such high energies, we also find that invariably 
only one or two diatoms are formed with the remaining atoms 
flying apart. That is, collision-induced dissociation occurs 
as is clear from the sample trajectories plotted in figs. 31 
and 32. 

We tried to follow the molecular collision in 3D by 
making different types of plots. In the absence of animation 
facilities in our computer, we plotted 'stills' of the six 
atoms at different Intervals of time. The composite diagram 
showing the time evolution became messy. Therefore we decided 
to show in one frame the molecules at the 'beginning' and at 
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Bond distance (X) 



Time (1 x lO^^s) 


Fig. 29 Bond distance versus time plot for a typical 
trajectory showing reaction ( R12) 



Bond distance (A) 



Fig. 30 Same as fig. 29 for a reaction 
of the type ( R 9)» 


Bond distance (A) 



Fig. 31 Same as fig- 29 with only two 
diatomic molecules as the 
products . 



Bond distance (A) 



Fig. 32 Same as fig. 29 with one stable diatomic 
molecule as a product, 


/y 


the 'end' oC the collision. A perspective view is shown in 
figs. 33 and 34 for the above two types of collisions 

We have investigated the effect of Eip and v on o in a 
liinlted way. The summary presented in fig. 3 5 shows that for 
(0,0,0), the threshold energy 90 kcal .mol~^, to be 

compared with the classical of 69 kcal.mol ^ for the 6C 
pathway. It is clear that there is a substantial increase 
in a with increase in F^p. At a fixed E^, a increases with 
V, Considering the fact that = 6.31, 18.38 and 31.20 
kcal.mol”"^ for v - 0,1,2 respectively, it is clear that energy 
partitioning in the vibrational mode of the attacking diatom 
enhances . The fact that, as for atom-diatom systems, the 

R T-l^ 

increase in d with E,j, arises mainly from an increase in P 
in each b-stream while the vibrational enhancement arises 
significantly from an increase in the b^^^^ is illustrated with 
the aid of (b) and 2iTb.P^(b) plots in fig. 36-39. 

We have computed the state-selected rate constant for 
the reactions (R9)-(R12) at 1000 K using the relation 

k(T) = <v.u'^> (^3) 

We interpolated a^(E;j,) graphically and obtained k as 

4.5xl0~^^ cm^ .molecule”^ s \ which is 3.0x10 times 

29 

smaller than the observed value of Bauer and Ossa for the 




Fig. 33 3D perspective view of reaction (Rll), 
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b (A) 


Fig, 36 Reaction probability as a function of 
impact parameter for different 
(OjO^v) at E-j-s 100 kcaLmofl 



I n b. H \ A 



b(A) 


Fig. 37 Partial cross section as a function of. 

b for different (0,0,v) at EjslOOkcal.moT. 




Fig. 38 Same as fig. 36 for v= 0,0,0 at 
four different Ej . 




2Trb.P 



Fig. 39 Same as fig. 37 for v = 0,0,0 at 
four different Ej . 




react Jon (R3) in the presence ot argon Despite the tact 
that a 6C pathway is energetically favored over the 4C path 
way, the tormer has a dynamical threshold of >90 Iccal.niol"^ 
nearly twice the experimental E of 42 kcal.mol”^. Thus 
our computed k is bound to be several orders of magnitude 

less than the experimental k. As a matter of fact, if an 

-1 -1 
oC 90 kcal.moL is used instead of 42 kcal.mol in the 

rate expression of Bauer and Ossa, k becomes comparable to 

ours emphasizing the role of the third body in the 4C 

exchange reaction. 



CHAPTER FIVE 


SUMMARY AND CONCLUSION 

Wg have tested the utilaty of various different analytic 

functions in fitting repulsive potential-energy curves for 
3 I- 

^2 ^ ^u^ ' NaLi. A comparison of the different 

analytacally fitted potentials for the test case of H 2 was 
presented in Chapter 2 with the help of 'In V vs r' and 
residual plots. We have extended the study to 3 -dimensions 
by fitting the limited amount of data available on the PES 
for the linear and bent geometries of H* (2pE’)in the 
configuration by the rotated-anti-Morse-curve-spline (RAMCS) 
approach. The quality of the fit was illustrated with the 
help of contour plots, for different apex angles. We conclude 
that an exponential -poly function gives the best fit of the 
repulsive potential for H 2 ( ^ud HeH and close to the best 

Cor 1162 and NaLi The traditional anti-Morse function does 
serve as a reasonable alternative, particularly in view of the 
fact that it can be generalized to two and three dimensions as 
is the case for (2pE' ) . 

We have tested the utility of the SM function in fitting 
the chemically accurate ab initio data of Siegbahn and Liu for 
the ground state PES for H^. Details of the fit along with a 
discussion of the quality of the fit with the aid of contour 
maps are given in Chapter 3. The ^ function has been shown to 



_ 1 

reproduce the surface v/xtb an rms deviation of 0.9 koal mol . 

A TDOM calculation for H + H 2 (v = 1) at 1000 K on the SM 
and SLTH surfaces vindicates the utility of the former in 
spite of the deficiencies like the tendency to show spurious 
wells in the region of the PES where there is a lack of ab 
initio data. 

We have also reported the analytic fit of the ab initio 

data of Naqer and Junqen for H^(2pE'), The SM function 

reproduces this surface with an rms deviation of 2,6 kcal. 

1 

mol . Unfortunately our atbempts to make the two surfaces 
degenerate for geometries have not been successful. We are 

pursuing this aspect further. 

Wo have presented the results of a dynamical study of 
a model six -center exchange reaction in Chapter 4. We have 
Loiind Ihdt at low (100-135 kcal.mol”^) the reaction pro- 
ceeds essent jolly as a 4C pathway with one of the diatomic 
mo3 ecules in the dimer acting as a mediating partner. At a 
higher E,p (13 5-155 kcal.mol"*^ and for large v = 8-11) a 'true' 
6C reaction takes place invariably leading to some dissociated 
molecules because of the high energies Involved The computed 
values yield a threshold energy of ~90 kcal. mol . While 
increases with E;p there is a vibrational enhancement when 
the energy is repartitioned from translation to vibration. We 
have computed the state— selected rate constant at 1000 K to be 
4.5x10“^^ cm^ molecule"^ several orders of magnitude 



smal Ler than the experi'tientaH y observed value This 
discrepancy arises largely fro'-n the fact that despite the 
6C pathway having a lower Ej^ than the 4C analog, the observed 
ol 90 heal. mol ^ is still a Cactor of two larger than 
the experimental E of ~ 42 hcal.mol”^. Since ours is the 

a. 

first 6C dynainical study, a much more elaborate study is 
needed and is planned for the next few months 
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Fitting Repulsive Potential-energy Curves and Sui faces 

ii\ SlJKAKMA IiIARIJAANdN Sai UVAMUIUUY’'' 

Dcptirimcnt ol Chcmisiiy, Inclmn Insltlulc onecimology, Kanpur 208 016, India 

Reiewed 29lh August, 1984 


The utility of various dilTorcnt analytic functionb in Pilling repulsive potential-energy 
curves lor a Pew diatomic species has been tested An exponential Punction with a polynomial 
exponent is louiui to be the best A limited nuinbei oP data on the potential-energy surfaces 
for the liito.ir aiul bent geoiiieiucs of Hf m the /2„, conhguralion) can be fitted 
leasomibly well by an appionch involving a rotated anti-Morse curve-spline interpolation 


hitimg i») oiiiilytic lunction to potential-energy values continues to remain a 
bottleneck in going Piom clcclionic-stnicture calculations to detailed dynamical 
atuclies. r-itting potentinl-cneigy vnlucs as a function of one independent variable 
IS the simplest pioblem. While consideiable progress has been made in fitting 
potential-energy cuives containing one minimum,’ the problem of fitting a purely 
repulsive ciiivc lias leccived less attention The fust attempt to represent a repulsive 
poiential-eneigy curve was by Sato^ and was for the tiiplet state C^u) oPHj He 
had used nn unti-Morse function of the foim 

V-»“[exp(“2/3p)-l-2cxp(-/3/j)] (I) 

where p =• r -le mid O, >c and (i arc the bond dissociation energy, the cquilibmim 
bond distance and the cuivauirc puiametcr, tespectively, for the giound state ('i.g} 
of lij and have no physical meaning foi the liiplet state. 

Pedersen and Poitei^ extended the function by matching it with an exponential 
function at longci langc* 

V « A[ exp ( ~2fip) -I 2 exp { -(ip)l r « 1 .6 au 

-l>(r-l-c)cxp(-ttr), r>l.6au (2) 

and were able to fit the ah initio potentials for II 2 and Kuntz and Roach® 

modified the anti-Morse function further by allowing /3 to be a variable. 

/3>=/3{l l-ypd 5p^) (3) 

and fitted the potentials for llj and II 2 successfully,*’’’ Berces“ has shown that 

V = (D/«)cxp(-2^p)[l + (/3p)’“] (4) 

IS also more flexible than the anti-Morsc function in fitting the triplet potential-energy 
curves of Hj, HI and I 2 . Kafrt and Berry* used a simple exponential with a 
polynomial exponent for fitting the ah Muo data of Kolos and Wolniewicz*' 

V » exp [(«o+ r + 02 ^+ <i 3 r’+ a 4 r‘'+ osr*)]. 
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REPULSIVE POTENTIAL ENERGY CURVES 


Stine and Muckerman’® used a one-dimensional spline*' for fitting the potential- 
energy cuive for Soniileitnei and Beckel'^ have shown that rational 

fractions can be used to reproduce the potential for H|(2/70’„) ” 

More recently Vaiaudas and Brandao'*' have proposed the useol an exponetitial 
multiplied by r“', i c 

V = /It exp ("h/j) (<i) 

for alkali-metal dimers in their '2. state Almost a decade ago Murrell and Sorbie'® 
had proposed the use of an extended Rydbeig function 

V = -D exp (-fl|p)(l + aip-fajp^+ ) (7) 

for fitting potentials with a single minimum However, there is no reason why it 
cannot be used to fit a repulsive curve m the same way as an anti-Moise iunction 
In principle, given adequate data, the potential-energy values can be fitted by a 
polynomial, or better still by a one-dimensional spline or Akima'^ interpolation, In 
practice one prefers a simple analytic function with a few parameters to be deter- 
mined 

To the best of our knowledge no comparative study has been made of the validity 
of these various lunctions in filling the data for a single system The obvious choice 
for a test case would be the stale of Ha for which accurate ab milio values are 
available"* We have undertaken such a study for Ha, HCj, HeH and NaLi and the 
results are summarized in the next section We must mention that we have only 
considered global functions in this study and not the use of piecewise interpolation 
methods We have specifically considered only the anti-Morse,* Berces,® exponen- 
tial-poly,** modified exponential '■* and extended Rydberg‘S functions, as we are 
interested in functions which require only limited input data and also are extendable 
to higher dimensions 

As a test case of a repulsive potential-energy surface with 2 or 3 independent 
variables we have examined the possibility ol fitting the potential-energy surface 
for H* (2pE' in its D 31 , geometry). This surface is needed in our study of the 
spectroscopy of the transition state m predicting the wings to the Lyman-a Ime in 
H + H2 collisions *’ The amount of ab tnuto information available on this surface 
is scanty *** Therefore, at our request Raynor'®* has computed a diatomics-in- 
molecules'®*' potential-energy surface in rotated Morse coordinates*® (/, 0) for linear 
as well as non-lincar geometries, with three potential-energy values along each 
diieciion, for d==0, 10, 20, 30, 35, 40, 42, 43, 44 and 45“, and hence by symmeliy 
at 0 = 46, 47, 48, 50, 55, 60, 70, 80 and 90“ for the bending angle y = 60, 90, 120 
and 180“. In the next section we report our successful rotated anli-Morse-curve- 
splinc (RAMCS) fit of the potential for fixed-y geometries 

RESULTS AND DISCUSSION 

Ab tntito values* of in V for H2(*X;I') are plotted against r in fig I, and for 
comparison the diflerent analytically fitted potentials are shown in the same graph. 
The parameters for each function were determined by a non-linear least-squares fit 
of the ab tnino data using the E04FCF program of NAO.*® The root-mean-squure 
(r.m s ) and the maximum deviations for each function are listed in table I, 

The anti-Morse function fits the H2 potential adequately at short range. The 
modification introduced by Berces improves the fit dramatically, but at the expense 
of adding a parameter The exponential-poly function deviates much more than 
the Berces function at long range However, the former gives the best overall fit, 
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Fig. 1. Compdnson of dilTereiU analytic fits to the afolHiRO data The aft mifio potential energy 
values' are connected by a solid line The dashed lines represent the different functional fits 
(1) modified exponential, (2) anti Morse, (3) extended Rydberg, (4) Bertes and (5) exponen 

tial-poly 


as IS illustrated with the aid of a residual plot in fig 2 The extended Rydberg 
functional hi is only slightly inferior to the Berces, but involves as many parameters 
as the exponcntial-poly fit The modified exponential function gives the poorest fit 
We have tried the different functional (us for the ab inido potentials for Hca, 
Hell and NaLi The i cs nils sunimarized in table I show that while the anti-Morse 

function fils the potential leasonably, the exponenlial-poly gives a better fit for all 
the systems The quality ol the In obtained by using the extended Rydberg function 
IS compaiable to that ol the exponential-poly for Hej and NaLi, Berces gives the 
woisl fit for NaLi and the extended Rydberg that for HeH 

If 11 large number ol potential-energy values are available and the number ol 
parameieis can be as large as 6, (he exponential-poly lit is obviously the best or 
close to being the best II also offers the advantage that In V can be fitted to a 
polynomial by a linear least-squares fit to yield an excellent set of initml-guess 
piiiamclcrs loi cairyiiig out a non-hneur Icasl-sqiiiues IK (an iinpoitaiil priKiieuI 
advantage over the other functional approaches) On the other hand, if the number 
ol pararneiets has to be as low as ihice, the anli-Morse (unction yields letisonable 
results and U also promises to be the candidate when potential-energy values 
involving more than one independeiu vaiiable are to be fiued 
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Fig. 2. Residual plots comparing the performance of the Berces (O), extended Rydberg (A) 

and exponential-poly (•) functions 

For example, for the collineui geometry of we have fitted the potential 
along each cliiection lor cliflerent values of 0 using an anti-Morse lunction The 0 
dependence ot the resulting parameters as plotted in hg 3 was reproduced by a 
one-dimcnsional spline interpolation method, and the overall RAMCS lit is illus- 
tiatcd with the aid of potential-energy contours in fig 4 We were also able to 
obtain an RAMCS fit to the potential for the bent geometries The d-dependcnce 
of the paiameieis for y « 120, 90 and 60“ is comparable to that of y = 180“ except 
lor a double hump in 1,(0) for y = 60“. The potenlial m the range 0-6 7 eV is 
reproduced with i.tns deviations of 0,21, 0 19, 0 15 and 020eY foi y— 180, 120, 
90 and 60“, respectively The maximum deviations are correspondingly 0 35, 032, 
0 34 and 0 30 eV When necessary, the ((?, y) dependence ol the anti-Morse para- 
meters cun be fitted by a two-dimensional spline" to obtain an overall three- 
dimensional fit 

We weie constrained to use only three-parameter functions in our routed- 
function approach, as only three potential-energy values were available along each 
direction, We tried to use 

Ysexp + 


and 


y-D/sxp (-a,0(fli + «20 


( 9 ) 
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Fig. 2. Residual plots comparing the performance of the Berces (O), extended Rydberg (A) 

and exponential-poly (•) functions. 


For example, for the collinear geometry of Hf we have fitted the potential'**^ 
along each direction for different values of 0 using an anti-Morse function. The d 
dependence of the resulting parameters as plotted in fig. 3 was reproduced by a 
one-dimensional spline interpolation method, and the overall RAMCS fit is illus- 
trated with the aid of potential-energy contours in fig. 4. We were also able to 
obtain an RAMCS fit to the potential for the bent geometries. The ^-dependence 
of the parameters for y = 120, 90 and 60“ is comparable to that of y = 180° except 
for a double hump in 4(0) for y = 60“. The potential in the range 0-6,7 eV is 
reproduced with r.m.s. deviations of 0,21, 0.19, 0.15 and 0.20 eV for y= 180, 120, 
90 and 60°, respectively. The maximum deviations are correspondingly 0.35, 0.32, 
0.34 and 0.30 eV. When necessary, the (0, y) dependence of the anti-Morse para- 
meters can be fitted by a two-dimensional spline" to obtain an overall three- 
dimensional fit. 

We were constrained to use only three-parameter functions in our rotated- 
function approach, as only three potential-energy values were available along each 
direction. We tried to use 


V = exp (ao+^iiP + <^ 2 P^) 

(8) 

V = £>/ exp (-ai/)(ai + a 2 /) 

(9) 
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which tire the truncated and modilieci forms of the exponential^poly and extended 
Rydberg functions, respectively We did not obtain satisfactory fits at low energies 
using eqn (8) Although cqn (9) gave a good ht along eacn direction, the variation 
of the pararneiers {D^ and « 2 ) with Q was not smooth, Hence we were not in a 
position to fit the potential successfully using these two functions 


CONCLUSIONS 

We have shown that, while an exponantial-poly function gives the best fit of the 
lepulsive potential lor and HeH, and close lo the best for Hca and NaLi, 

the traditional anti-Morse function does serve as a reasonable ahernaiive, par- 
ticuUuly m view of the lacl that it can be generalized to two and three dimensions 
as for {IpE^) 

Note udded w proof Huxley et have reported recently that an appended 
extended Rydberg liinclion 

4(1 + -f exp(-a,/ ) 

-U\nh(i -i„,/2)(c,,r ^*'h<8r i c.^r 

performs beltei than the extended Rydbeig function in luting repulsive potentials with 
n van der Waa\s minimum lor a Luge niimbci of diaioniic species 

The calculations reporlcd hcic were cairied out on a DCC-I()9() Computer at 
the Indian Institute ol Technology, Kanpur I his study was suppoited in part by 
a grant liom the Depailmenl ol StiencLand lechnology, New Dclin Weait tliankliil 
lo the lekiees loi tlieii comments 
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